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We study the case of "PT-symmetric perturbations of Hermitian Hamiltonians with degenerate 
eigenvalues using the example of a triple-well system. The degeneracy complicates the question, 
whether or not a stationary current through such a system can be established, i.e. whether or not 
the PT-symmetric states are stable. It is shown that this is only the case for perturbations that 
do not couple to any of the degenerate states. The physical explanation for the inhibition of stable 
currents is discussed. However, introducing an on-site interaction restores the capability to support 
stable currents. 
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I. INTRODUCTION 

Bose-Einstein condensates (BECs) of dilute ultracold 
atom gases are one of the most prominent examples of 
many-particle systems that can be described by a mean- 
field theory. For temperatures that are considerably 
smaller than the critical temperature of the condensate 
the Gross-Pitaevskii equation [1, 2] is known to provide 
very accurate results, although the description has limita¬ 
tions in the vicinity of dynamic instabilities [3, 4] . 

There is a broad field of applications of such BEGs in 
fundamental research, covering anything from the sim¬ 
ulation of condensed matter in optical lattices over the 
study of rotating ultracold gases to the transition from 
fermionic superfluids to BECs [5] . Another promising ap¬ 
plication is the atom laser [6]. To realize a stable source 
of coherent atoms for the study of transport phenomena 
or the realization of an atom laser, particles have to be 
coupled in and out of the condensate coherently. Both 
has already been realized experimentally [7, 8]. 

To make use of the Gross-Pitaevskii equation even in the 
case of particle exchange with the environment, one can 
turn to non-Hermitian systems [9]. Here particle in- and 
outflux are described by positive and negative imaginary 
potentials, respectively [10]. Exact many-particle calcula¬ 
tions [11] and master equation approaches [11-14] support 
the use of the non-Hermitian mean-field ansatz. Other 
applications of non-Hermitian Hamiltonians in dissipative 
optical lattices [15, 16] and in the study of condensate de¬ 
cay [17, 18], as a result of a complex scaling approach [19], 
have also been very successful. 

However, it is not quite clear whether these Hamiltoni¬ 
ans provide really stationary states supporting a stable 
current through the system since the energy eigenvalues 
can take complex values. A condition that is known to 
provide real eigenvalues in many cases is TT symme¬ 
try [20-22]. Such systems have either real eigenvalues, 
belonging to TT-symmetric eigenstates, or complex con¬ 
jugate eigenvalues, belonging to TT-broken eigenstates. 
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The concept of VT symmetry is a special case of the a the¬ 
ory of pseudo-Hermitian operators [23-25] , however, both 
properties become identical in the case of a Schrodinger 
equation in position space with complex potentials. 

The capability of these systems to provide real eigen¬ 
value spectra was analyzed in many theoretical studies [26- 
33] and was realized experimentally in optical wave-guide 
systems [26, 34-36], which promises technical applica¬ 
tions. Even though proposals exist [37, 38], the concept 
has not yet been realized experimentally in ultra-cold 
atomic gases. 

These concepts have in common that a Hermitian 
system is perturbed by a non-Hermitian TT-symmetric 
Hamiltonian. For weak perturbations such systems can 
be discussed using perturbation theory. This technique 
has already been used with great success [39, 40] to de¬ 
cide whether a linear system stays VT symmetric under 
a non-Hermitian perturbation. Since complex eigenvalues 
lead to an exponential growth or decay of the eigenstate, 
TT-symmetric states in a linear system are stable if and 
only if the eigenvalue spectrum is entirely real, i.e. the 
symmetry is unbroken. 

We want to tackle the question under which conditions 
a system, in particular a three-well setup, is capable of 
supporting a stable stationary current, i.e. stable station¬ 
ary states for finite TT-symmetric imaginary potentials. 
For simple, i.e. non-degenerate, discrete eigenvalues which 
become equally spaced for high quantum numbers, it has 
already been shown that real eigenvalues always exist and 
can be found if the non-Hermitian part of the Hamilto¬ 
nian is small enough [41-43]. However, this situation can 
change completely if degenerate eigenvalues are considered. 
The focus of the present work lies on the question how 
such degeneracies influence the existence of stable VT- 
symmetric states. Additionally, since we are interested in 
BECs, the impact of the nonlinear contact interaction in 
the Gross-Pitaevskii equation on the linear results will be 
studied. In this case the stability has to be investigated 
separately using the Bogoliubov-de Gennes equations. 

In Sec. H, we recapitulate the perturbation theory ap¬ 
proach for both anti-Hermitian and pseudo-Hermitian per¬ 
turbations and derive a stringent condition under which 
the VT symmetry is preserved at least for smallest per- 
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turbations. Then we turn to a three-well system whose 
properties are discussed in the linear and nonlinear case 
in Sec. Ill, testing the prediction from the perturbation 
theory approach and the influence of the inter-particle 
interaction. Finally conclusions are drawn in Sec. IV. 

II. PERTURBATION THEORY 

We discuss a linear system described by the Hamilto¬ 
nian H = Hq + 'yHp , where the unperturbed Hamilto¬ 
nian Hq is Hermitian and the perturbation Hp is anti- 
Hermitian. If they fulfill the relations 

VHo = hIv = HoV, (la) 

VHp = hIv = -HpV, (lb) 

the Hamiltonian is pseudo-Hermitian with respect to the 
parity operator, VH — WV. It is obvious that in case of 
a diagonal matrix representation, the system is VT sym¬ 
metric, where T acts as a complex conjugation. This is in 
particular true in position space where the perturbation 
is given by an antisymmetric imaginary potential iV (x). 

If all eigenvalues of Hq are non-degenerate and discrete 
a perturbation can always be chosen small enough such 
that the eigenvalues, provided they behave as continuous 
functions of the perturbation strength 7 , do not coalesce 
and form exceptional points. Since complex eigenvalues of 
pseudo-Hermitian Hamiltonians always appear in complex 
conjugate pairs [23-25] the formation of complex eigenval¬ 
ues can only occur at exceptional points. In particular, for 
perturbations small enough all eigenvalues remain real. 

This particular property can also be seen in a more 
formal way using perturbation theory [44, 45[. In ad¬ 
dition, this approach can later be used to analyze de¬ 
generate systems. Since the unperturbed Hamiltonian 
commutes with the parity operator V, all eigenvectors, 
{IV'm)}, can be chosen either symmetric or antisymmet¬ 
ric, hence Vlipm) = = ±1 holds. The power 

series for the energy eigenvalue of the corresponding 
eigenstate ji/'n) reads 

^J■n = Mn.o + 7 /^ 11,1 + 7^Mn,2 

A useful explicit form of the energy corrections was given 
by Tosio Kato in 1949 [46[, 

= (- 1 )"-^ tr ^ S'^^HpS^^ ■ ■ ■ HpS ’^‘+^, ( 2 ) 

ks 

where the sum is chosen to iterate over all G with 
A:i + ... + ks-\-i = s — 1. The reduced resolvent is defined 
as 

IJ-m- fj-n’ 

m^n 

while the zeroth power is treated differently as 
S'" = \n){n\. 


We see that [S, V] = 0 since 

7^|m)(m| = Xm.\m){m\ = \m){m\Xm = |m)(m|7^. 

Pseudo-Hermiticity of Hp ensures that all corrections 
fj,n,s are real. This can be seen by using tr{VOV) = tr(0) 
and = S and calculating the complex conjugate of the 
energy correction 

= {-ly-HiY.rs'^^+^Hl ■ ■ ■ hIs'^^v 

ks 

= (-l)®-4r^S'==+iiSpP---PHpS'=i 

ks 

= (-l)®-4r^S'=iHp---iSpS'==+i 

ks 

— (3) 

However, since anti-Hermiticity of Hp renders every odd 
energy correction purely imaginary, 

rk,s = (-!)*■' tr ^ ... HlS^^ 

k^ 

= (- 1 )®”^ H^i-iys’^^Hp ■ ■ ■ HpS^‘+^ 

k^ 

= (-l)>n.., (4) 

they vanish and we conclude that the lowest order energy 
shift is at least quadratic. 

The situation changes drastically if we stop assuming 
non-degenerate eigenvalues. Let us consider M degenerate 
eigenvectors {4>i} with the common eigenvalue /i. The 
first order energy shift is given by the eigenvalues of the 
matrix 

Sij = {(t)^\Hp\(|)J). (5) 

The matrix S inherits the anti-Hermiticity from the per¬ 
turbation Hp. Its eigenvalues are therefore either purely 
imaginary or zero. This means that any perturbation that 
does not imply S' = 0 immediately produces imaginary 
eigenvalues. For PT-symmetric Hamiltonians describing 
in- and outcoupling of particles, this means that no stable 
stationary current from particle source to particle sink 
can exist. 

For S = 0, i.e. all matrix elements between degenerate 
eigenstates are zero, higher-order perturbations must be 
taken into account. In this case the considerations for 
simple eigenvalues from Eq. (3) hold, thus all eigenvalues 
remain real. 


III. THREE-MODE SYSTEM 

To gain a more physical understanding of the effects 
leading to this symmetry breaking that prohibits any 
stable stationary current through a system, we turn to 
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Figure 1. Scheme of the three-mode model. The imaginary 
on-site energies provide a particle influx to well 1 and an 
outflux from well 3. These two wells are coupled by a constant 
strength, set to unity, while the coupling J to well 2 is varied. 
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a simple example and solve the three-mode GPE with 
on-site interactions, 

3 

+ U \ipi\'^'ipi = ( 6 ) 

i=i 

described by the Hamiltonian 

/ 0 -J -l\ /i 0 0 \ 

H = - J 0 - J -b 7 0 0 0 . (7) 

\-l -J 0 J \0 0 -ij 

'■ -V-" '-V-" 

Ho Hp 

The parameter U is the strength of the nonlinear on-site 
interaction. The chemical potential /r takes the role of 
the energy and, for the linear case, C/ = 0 , is given by 
the eigenvalues of the Hamiltonian. Therefore the three 
terms are used as synonyms. 

The Hamiltonian describes a system of three coupled 
wells where in the first well particles are injected into the 
system while particles are removed from the third well 
at the same rate 7 . The tunneling rate between the gain 
and loss well is fixed to unity while the second well is 
coupled to both by the same variable coupling strength 
J. A scheme of the setup is shown in Fig. 1. 

It is obvious that this Hamiltonian is PT-symmetric 
and pseudo-Hermitian with respect to 

/O 0 1 \ 

P = 0 1 0 . 

Vl 0 oj 

The effects leading to the symmetry breaking intro¬ 
duced in the previous section are discussed in the linear 
system, U = 0. Without the perturbation, i.e. 7 = 0, the 
eigenvalues are given by the simple expressions 

/ri = -V 2 J 2 + l/4-l/2, ( 8 ) 

/r 2 = +x/2J2 + l/4-l/2, (9) 

= 1 - ( 10 ) 

There are two cases that deserve special attention. For 
J = 0 the second well is not coupled to the remaining 


Figure 2. (Color online) (a) Real and (b) imaginary part of the 
eigenvalue /r for different values of the coupling strength J as 
function of the in- and outcoupling parameter 7 . For J = 0 the 
ground and second excited state form an exceptional point of 
second order (EP2) at 7 = 1 where both "PT-symmetric states 
vanish and two "PT-broken states emerge. The first excited 
state, located only in well 2 is completely independent from 
the other eigenstates. For J > 0 the first excited state forms 
an EP2 with the second excited state, while the ground state 
is shifted to smaller energies and retains its PT symmetry for 
all values of 7 . For increasing values of J, the energy difference 
between the two excited states decreases and the EP2 is shifted 
to smaller values of 7 . At J = 1, two PT-broken states emerge 
directly at 7 = 0 with degenerate eigenvalues jj, = 1. Thus 
rendering the ground state unstable for 7 7 ^ 0 . 

double-well system. This case has already been discussed 
in detail [33]. The PT symmetry in the system breaks 
at 7 = 1. The second special case is J = 1, in which 
the system becomes totally symmetric. In this situation 
the two eigenvalues P 2 = M 3 = 1 degenerate. Applying 
perturbation theory sketched in the previous section to 
this case leads to the matrix 

»=(2io)' (“> 

immediately breaking the PT symmetry. 

The transition from the double well to the symmetric 
triangular case is now studied in more detail. The com¬ 
plete analytic expressions for 7 7 ^ 0 do not provide further 
insight, therefore we limit our discussion to numerical 
results starting with the linear spectrum in Fig. 2. For 
J = 0, where the system corresponds to a double-well 
potential and one independent well, the ground and ex¬ 
cited state of the double well coalesce in an exceptional 
point of second order (EP2) where they vanish and two 
PT-broken states emerge. These states have complex 
conjugate eigenvalues with a finite imaginary part. The 
remaining PT-symmetric state is completely localized 
in the decoupled well 2 , while the double well is empty. 
However, any perturbation that adds some number of 
particles to the wells 1 and 3 still renders the remaining 
double-well system unstable. 

For J > 0 the structure changes and the ground state 
stays PT-symmetric for all values of 7 . The EP2 however 
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Figure 3. (Color online) Net current through the system jext 
as function of 7 for (a) J — 0, (b) J = 0.1, (c) J = 0.4, (d) 
J = 0.8. The highlighted areas show regions where "PT-broken 
eigenstates with complex eigenvalues exist, thus rendering all 
stationary states unstable. The maximum current, jext = 1, is 
reachable only for J = 0. For 0 < J < 1 the second excited 
state supports the highest net current slightly before the EP. 


is still present, now occurring between the two excited 
states of the three-well system. For increasing coupling 
strengths J the EP2 is shifted to lower parameters 7 
breaking the VT symmetry of these states and rendering 
the system unstable at smaller values of 7 . 

For J = 1 the EP2 reaches 7 = 0 where it vanishes 
due to the Hermiticity of the system and is replaced by a 
degeneracy of two eigenvalues. For this case, as predicted 
from perturbation theory, the VT symmetry is broken 
for any value of 7 . 

The coupling to an additional well and thus the avail¬ 
ability of an additional channel seems to diminish the 
capability of the system to support a stable current. The 
net current through the system is given by the particle 
current from an external source to the gain well 1. The 
continuity equation for this well reads 


^ |ui|^ = 27 |ui|^ - Y {vlv 2 - V1V2) - Y - U1U3), 
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Jl3 


( 12 ) 

where Vi are the components of the eigenvector and with 
the particle currents jij from well i to j and the external 
current of particles into the well jext- For stationary states 
jext = J12 + J13 holds. At the same time the incoupling 
in well 1 and the outcoupling in well 3 are balanced, thus 
jext specifies the net current. 

Figure 3 shows the net flow jext for all three states and 
different coupling strengths J. For J = 0 the ground 
state and second excited state, which are the solutions to 
an equivalent two-mode problem, support equally strong 
currents. The maximum current jext = 1 is reached at 
the exceptional point at 7 = 1. If the additional well 2 
is coupled with 0 < J < 1, the ground state is no longer 
involved in the EP2 and provides a stationary current even 
after the two excited states have vanished. However in 
this region the state is no longer stable with respect to any 


Figure 4. (Color online) Ratio of the current 712 through the 
additional well and ji 3 directly to the loss well as function of 
7 for (a) J = 0, (b) J = 0.1, (c) J = 0.4, (d) J = 0.8. While 
the current ji 3 is always positive, negative ratios represent 
negative currents J 12 . Again the highlighted areas show regions 
where PT-broken eigenstates exist. For J = 0 no particles can 
be transferred through well 2. After coupling the well to the 
system the ground state supports a positive current through 
the well, while the excited states show negative currents. For 
J —> 1 the ratio approaches —1. 


small perturbations. The strongest current is achieved by 
the second excited state for a value of 7 slightly smaller 
than the EP2. The VT symmetry breaking at the EP2 is 
therefore the main reason for the decrease of the maximal 
currents as the coupling J is increased. 

As a next step we look at the partial particle flows 
from the gain well 1 to the loss well 3, jis, i.e. the direct 
channel, and to well 2 , ji 2 , i.e. the additional channel. 
These currents are compared in Fig. 4 for different cou¬ 
pling strengths J. The ground state is the only stationary 
state supporting a positive current through the additional 
channel. For increasing coupling parameters J an in¬ 
creasing part of the particles - for J = 0.8 (Fig. 4(d)) - 
approximately a third of the net current, is transported 
through well 2 . 

In contrast to this behavior the excited states do not at 
all support a positive current through well 2. Since the 
current J 12 is negative, the additional channel transports 
particles from the loss to the gain well. The relative 
strength of this reverse current increases with higher 
parameters J, and becomes comparable to the current 
from the gain to the loss well jis for J —)■ 1 (Fig. 4(c), (d)). 
For J « 1 a small parameter 7 is sufficient to induce a 
strong circular current, breaking the VT symmetry even 
though only very few particles enter the system. 

Up to now we only discussed linear systems, i.e. sys¬ 
tems without inter-particle interaction. However, in most 
Bose-Einstein condensates at least a contact interaction, 
described by the on-site interaction in Eq. ( 6 ) is present. 
In the nonlinear case, the stability of a stationary state 
■0 is no longer determined by the complex eigenspectrum 
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Figure 5. (Color online) Chemical potential /r of VT- 
symmetric states against the strength of the parameter 7 
for coupling parameter J = 1 and nonlinearity parameter 
(a),(b) 17 = 1, (c) ?7 = 2, (d) 17 = 4. Introducing the non¬ 
linearity (a) gives rise to four new states for 7 7 ^ 0 , which 
are shown in the enlarged figures (b),(c),(d). For 17 = 1 and 
U = 2 three of them are stable (solid lines) and one is unstable 
(dashed line), for 17 = 4 all states except the highest excited 
one are unstable. 

but by the linear Bogoliubov-de Gennes equations: 

3 

HijUj + 2U Ui + Uiljfvi — ^Ui = ujUi, (13a) 
i=i 
3 

H*jVj + 2U - ^*Vi = -ojVi, (13b) 

i=i 

defining the linearized evolution of the perturbation 
S-ipi = Mi exp (—iwt) -b M* exp (iw*<). If all eigenvalues 
oj are real, the state is considered to be stable, while any 
non-vanishing imaginary part renders it unstable. 

Figure 5 shows the chemical potential for J = 1 and 
different strengths of the nonlinearity U. Again only VT- 
symmetric states are shown. The first point to notice is 
that every nonlinearity strength gives rise to a set of four 
new T^T-symmetric states (Fig. 5(a), (b)). Since three of 
these states are stable, the particle interaction supports 
a stationary transport through the system although the 
symmetry of the system inhibits all stable currents in 
the linear case. The four states vanish in two indepen¬ 
dent tangent bifurcations. For a stronger nonlinearity 
U = 2 (Fig. 5(c)) two stable states are available up to 
7 = 1. For U = A (Fig. 5(d)) only the highest excited 
stationary state is still stable. As already known from in 
T^T-symmetric systems the stability of stationary states 
does not change exactly at a bifurcation point but only 
close to them. However, in the following discussion we 
neglect this subtlety. 

The three states supporting stationary currents are 
now discussed in more detail. For that purpose, Fig. 6 
shows the particular net current for different strengths 
of the nonlinearity. For lower nonlinearities the three 
stable states support a similar maximum net current 
but at different values of 7 (Fig. 6 (a)). This is due to 


Figure 6. (Color online) Net current through the system given 
by the incoupling external current jext against 7 for (a) 1 / = 1 , 
(b) U = 1.5, (c) U = 2, (d) U = 2.5. Only the three stable 
"PT-symmetric states introduced in Fig. 5 are shown. For 
nonlinearities of U up to 2, all three states are able to sustain 
a similar current, which reaches a maximum of approximately 
0.3 at 1/ « 2. For higher nonlinearities stable states become 
unstable until only the highest excited state remains. Its low 
current is due to the fact that there are only few particles 
present in the gain and loss well. 


different values of the modulus square of the different 
wave functions at wells 1 and 3. The second excited state 
is almost fully localized in these wells. This explains why 
this state vanishes for lower values of 7 but is already 
supporting a strong current. For U ~ 2 the currents reach 
their maximum of about 0.3, i.e. a third of the maximum 
current of the double well with J = 0. For higher values 
of U the first and second excited state become unstable. 
The fourth excited state however stays stable since the 
weak population of the gain and loss well leads to a weak 
net current. 


IV. CONCLUSION AND OUTLOOK 

Perturbation theory predicts that in linear systems a 
PT-symmetric perturbation i7p of a Hermitian system 
leads to real eigenvalues and therefore stable stationary 
states if and only if all sets of degenerate eigenvectors 
{(pi} satisfy {(p,\HY>\(t>j) = 0 . 

As a simple example we studied a triple well with one 
gain, one loss, and one additional well. This system 
violates the above condition in the case of equal coupling 
strengths between all wells. The physical nature of the 
symmetry breaking was discussed by varying the coupling 
strengths while preserving the VT symmetry. It was 
demonstrated that the current supported by the system 
reaches its maximum when the third well is not coupled 
to the system, ergo in the case of a double well system. 
Coupling the well to the system leads to a reverse current, 
reducing and ultimately quenching the net current. 

While the equally coupled triple well does not support 
stable currents in the linear case, introducing a contact 
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interaction gave rise to four new T^T-symmetric states, 
three of them stable. These states restore the capability 
of the system to support a stationary current. For an 
interaction strength above twice the coupling strength 
between the wells, the current reaches its maximum of a 
third of the double-well current in the linear system. 

It would be highly interesting for future studies to ana¬ 


lyze spatially extended three-dimensional triple wells, i.e. 
a full solution in the three dimensional position space. 
While the results from perturbation theory provide good 
predictions for linear systems, it is not clear if the inter¬ 
particle interaction in extended potentials is able to re¬ 
store the support for stationary currents. 
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